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PUSHED BEYOND THE BRINK: ALLEE EFFECTS, ENVIRONMENTAL 
STOCHASTICITY, AND EXTINCTION 



To appear in a special issue on Allee effects in the Journal of Biological Dynamics 

Abstract. To understand the interplay between environmental stochasticity and Allee effects, 
we analyze persistence, asymptotic extinction, and conditional persistence for stochastic difference 
equations. Our analysis reveals that persistence requires that the geometric mean of fitness at low 
densities is greater than one. When this geometric mean is less than one, asymptotic extinction 
occurs with high probability for low initial population densities. Additionally, if the population 
only experiences positive density-dependent feedbacks, conditional persistence occurs provided the 
geometric mean of fitness at high population densities is greater than one. However, if the population 
experiences both positive and negative density-dependent feedbacks, conditional persistence only 
occurs if environmental fluctuations are sufficiently small. We illustrate counter-intuitively that 
environmental fluctuations can increase the probability of persistence when populations are initially 
at low densities, and can cause asymptotic extinction of populations experiencing intermediate 
predation rates despite conditional persistence occurring at higher predation rates. 



Populations exhibit an Allee effect when at low densities individual fitness increases with den- 
sity [Allee, 1931, Stephens et al., 1999]. Common causes of this positive density-dependent feedback 
include predator saturation, cooperative predation, increased availability of mates, and conspecific 
enhancement of reproduction [Courchamp et al., 1999, Stephens et al., 1999, Gascoigne and Lip- 
cius, 2004, Courchamp et al., 2008, Gascoigne et al., 2009, Kramer et al., 2009]. When an Allee 
effect is sufficiently strong, it can result in a critical density below which a population is driven 
rapidly to extinction through this positive feedback. Consequently, the importance of the Allee 
effect has been widely recognized for conservation of at risk populations [Dennis, 1989, Stephens 
and Sutherland, 1999, Berec et al., 2007, Courchamp et al., 2008] and management of invasive 
species [Keitt et al., 2001, Leung et al., 2004, Tobin et al., 2011]. Population experiencing envi- 
ronmental stochasticity and a strong Allee effect are widely believed to be especially vulnerable 
to extinction as the fluctuations may drive their densities below the critical threshold [Courchamp 
et al., 1999, Dennis, 2002, Berec et al., 2007, Courchamp et al., 2008]. However, unlike the deter- 
ministic case [Cushing, 1988, Gyllenberg et al., 1996, Schreiber, 2003, Yakubu, 2003, Jang, 2006, 
Luis et al., 2010, Elaydi and Sacker, 2010, Kang and Lanchier, 2011, Duarte et al., 2012, Kang and 
Udiani, 2014], the mathematical theory for populations simultaneously experiencing an Allee effect 
and environmental stochasticity is woefully underdeveloped (see, however, Dennis [2002]). 

To better understand the interplay between Allee effects and environmental stochasticity, we 
examine stochastic, single species models of the form 



where X t £ [0, oo) is the density of the population at time t, / (x, £) is the fitness of the population 
as a function of its density and the environmental state £, and the environmental fluctuations £t are 
given by a sequence of independent and identically distributed (i.i.d.) random variables. Here we de- 
termine when these deterministic and stochastic forces result in unconditional stochastic persistence 
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(i.e. the population tends to stay away from extinction for all positive initial conditions with prob- 
ability one), unconditional extinction (i.e. the population tends asymptotically to extinction with 
probability one for all initial conditions), and conditional stochastic persistence (i.e. the population 
persists with positive probability for some initial conditions and goes extinct with positive probabil- 
ity for some, possibly the same, initial conditions). Section 2 describes our standing assumptions. 
Section 3 examines separately how negative-density dependence and positive-density dependence 
interact with environmental stochasticity to determine these different outcomes. For models with 
negative density-dependence (i.e. is a decreasing function of density x), [Schreiber, 2012] 

proved that generically, these models only can exhibit unconditional persistence or unconditional 
extinction. For models with only positive density-dependence (i.e. f(x, £) is an increasing function 
of density x), we prove that all three dynamics (unconditional persistence, unconditional extinction, 
and conditional persistence) are possible and provide sufficient and necessary conditions for these 
outcomes. Section 4 examines the combined effects of negative and positive density dependence on 
these stochastic models. We prove that conditional persistence only occurs when the environmental 
noise is "sufficiently" small. Throughout all of the sections, we illustrate the main results using 
models for mate limitation and predator saturation. Section 5 concludes with a discussion of the 
implications of our results, how these results relate to prior results, and future challenges. 

2. Models, assumptions, and definitions 

Throughout this paper, we study stochastic difference equations of the form given by equation (1). 
For these equations, we make two standing assumptions 

Uncorrelated environmental fluctuations: {^}^ 0 is a sequence of independent and iden- 
tically distributed (i.i.d) random variables taking values in a separable metric space E (such 



Fitness depends continuously on population and environmental state: the fitness func- 
tion / : M + x E — > is continuous on the product of the non-negative half line M+ = [0, oo) 
and the environmental state space E. 

The first assumption implies that (Xt)t>o is a Markov chain on the population state space IR+. 
While we suspect our results hold true without this assumption, the method of proof becomes 
more difficult and will be considered elsewhere. The second assumption holds for most population 
models. 

Our analysis examines conditions for asymptotic extinction (i.e. limj_ s>00 Xt = 0) occurring with 
positive probability and persistence (a tendency for populations to stay away from extinction) with 
positive probability. Several of our results make use of the empirical measures for the Markov chain 
(X t )t> 0 given by 



where S x denotes a Dirac measure at the point x i.e. 5 X (A) = 1 if x £ A and 0 otherwise. For any 
interval [a, b] of population densities, lit ([a, b]) is the fraction of time that the population spends in 
this interval until time t. The long-term frequency that (Xt)t>o enters the interval [a, b] is given by 
lim^oo IIf([a, &]), provided the limit exists. As these empirical measures depend on the stochastic 
trajectory, they are random probability measures. 



3.1. Results for negative-density dependence. For models with only the negative density 
dependence (i.e. fitness / is a decreasing function of density), Schreiber [2012] proved that the 
dynamics of the model (1) exhibit one of three possible behaviors: asymptotic extinction with 





s=0 



3. Negative- versus positive-density dependence 
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probability one, unbounded population growth with probability one, or stochastic persistence and 
boundedness with probability one. Closely related results have been proven by Chesson [1982], 
Ellner [1984], Gyllenberg et al. [1994], Fagerholm and Hognas [2002] and Vellekoop and Hognas 
[1997]. Prior to stating this result, recall that log + x = maxjloga;, 0}. 

Theorem 3.1 (Schreiber 2012). Assume f(x, £) is a positive decreasing function in x for all £ £ E 
and E[log+ /(0, &)] < oo. Then 

Extinction: j/E[log /(0, < 0, then linn._ s . 00 Xt = 0 with probability whenever Xq = x > 0, 
Unbounded growth: if linix-^oo E[log f(x, £t)] > 0, then limt_> 00 = oo with probability 

whenever Xq = x > 0, and 
Stochastic persistence: if E[log /(0, £)] > 0 and ]im x -^oo E[log f(x, £t)] < 0, then for all 

e > 0 there exists M > 0 such that 

Iimsupllt([l/M, M]) > 1 — e almost surely 

t— >oo 

whenever Xq = x > 0. 

In the case of stochastic persistence, Theorem 3.1 implies that the typical trajectory spends most 
of its time in a sufficiently large compact interval excluding the extinction state 0. 

To illustrate Theorem 3.1, we apply it to stochastic versions of the Ricker and Beverton-Holt 
models. For the stochastic Ricker model, the fitness function is f(x, £) = exp(r — ax) where 
£ = (r, a). Stochasticity in r t and at may be achieved by allowing r t to be a sequence of i.i.d. 
normal random variables or at to be a sequence of i.i.d. log-normal random variables. These 
choices satisfy the assumption linx^oo E[log f(x, £t)] = — oo. This stochastic Ricker model is almost 
surely persistent if E[log/(0, £t)] = E[r^] > 0. If E[r^] < 0, then asymptotic extinction occurs with 
probability one. 

For a stochastic version of the Beverton-Holt model, we have f(x,£) = a/(l + bx) with £ = 
(a,b). Stochasticity in at and b t may be achieved by allowing them to be sequences of i.i.d. log- 
normal random variables. These choices satisfy the assumption hin^oo E[log f(x, £t)] = — oo. This 
stochastic Beverton-Holt model is almost surely persistent if E[log /(0, &)] = E[logdf] > 0. If 
E[loga t ] < 0, then asymptotic extinction occurs with probability one. 

3.2. Results for positive density-dependence. In contrast to models with only negative density- 
dependence, models with only positive density-dependence exhibit a different trichotomy of dynam- 
ical behaviors: asymptotic extinction for all initial conditions, unbounded population growth for 
all positive initial conditions, or conditional persistence in which there is a positive probability of 
the population going asymptotically extinct for some initial conditions and a positive probability 
of unbounded population growth for some, possibly the same, initial conditions. To characterize 
this trichotomy, we say {0, oo} is accessible from the set B C (0, oo) if for any M > 0, there exists 
7 > 0 such that 

P [{3t > 0 : Xt £ [0, 1/M] U [M, oo)} | X 0 = x] > 7 

for all x G B. 

Theorem 3.2. Assume f{x^) is an increasing function of x for all £ € E. Define foo(0 = 
lim^oo f(x,£). Then 

Extinction: if E[log /oo(^t)] < 0, then lim^oo Xt = 0 with probability one whenever Xq = 
x>0. 

Unbounded growth: if E[log /(0, £t)] > 0, then limt->oo X t = 00 with probability one when- 
ever Xq = x > 0. 
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Figure 1. Effect of initial population density on persistence for the stochastic mate limitation 
model. The stochastic mate limitation model with f(x, £) = where £ = (A, h) was simulated 
10, 000 times for each initial density. The fraction of runs where the final density was greater than 
100 are plotted as a function of initial density xo- Parameters: h = 10 and A log-normally 
distributed with log-normal mean 0.1 and log-normal standard deviations a as shown. 



Conditional persistence: i/E[log/(0, &)] < 0 and E[log/oo(£t)] > 0, then for any 0 < 5 < 
1, there exist m, M > 0 such that 



lim Xt = oo 



X n 



x 



> 1-5 and P 



lim X t 

t— >co 



Xn 



> 1-5, 



for all x G [M, oo) and a// y G (0, m]. 
Moreover, if {0, oo} is accessible, then 



lim X t 



0 or oo 



X 0 = x 



/or aZZ a; G (0, oo) 



To illustrate Theorem 3.2, we apply it to stochastic versions of models accounting for mate- 
limitation and a predator-saturation. For many sexually reproducing organisms, finding mates 
becomes more difficult at low densities. For instance, pollination of plants by animal vectors 
becomes less effective when patches become too small because lower densities result is reduced 
visitation rates by pollinators [Groom, 1998]. Alternatively, fertilization by free spawning gametes 
of benthic invertebrates can become insufficient at low densities [Knowlton, 1992, Levitan et al., 
1992]. To model mate limitation, let x be the density of females in the population. Assuming a 
50-50 sex ratio (i.e. x also equals the density of males in the population), Dennis [1989], McCarthy 
[1997], Scheuring [1999] modeled the probability of a female finding a mate by the function 



h + x 
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where h is a half-saturation constant i.e. the male density at which 50% of the females find a mate. 
If A is the number of daughters produced per mated female, then the fitness function is 

Ax 

0 = IT, — where £ = (A, h) 
h + x 

Stochasticity in £ t may be achieved by allowing At, ht to be sequences of i.i.d. log-normally dis- 
tributed random variables. Since E[log /(0, &)] = — oo, this stochastic model always exhibits as- 
ymptotic extinction for some initial conditions with positive probability. Theorem 3.2 implies that 
asymptotic extinction occurs for all initial conditions with probability one if E[lim x _ s>00 /(x,£t)] = 
E [log Xt] < 0. On the other hand, conditional persistence occurs if E [log At] > 0. 

Figure 1 illustrates how the probability of persistence for the mate limitation model depends 
on initial condition and the level of environmental stochasticity. Interestingly, higher levels of 
environmental stochasticity promote higher probabilities of persistence when initial population 
densities are low. Interestingly, when the population is below the "Allee threshold" , environmental 
stochasticity provides opportunities of escaping the extinction vortex. 

Another common Allee effect occurs in species subject to predation by a generalist predator 
with a saturating functional response. Within such populations, an individual's risk of predation 
decreases as the population's density increases. For example, in field studies, Crawley and Long 
[1995] found that per-capita rates of acorn loss of Quercus robur L. to invertebrate seed predators 
were greatest (as high as 90%) amongst low acorn crops and lower (as low as 30%) on large acorn 
crops. To model Allee effects due to predator saturation, Schreiber [2003] used the following fitness 
function 

f(x, £) = exp I r — I where £ = (r, P, h) 



h + x y 

where r is the intrinsic rate of growth of the focal population, P is the predation intensity, and h is 
a half-saturation constant. Stochasticity may be achieved by allowing rt to be normally distributed 
and ht,Pt be log-normally distributed. Theorem 3.2 implies that unbounded growth occurs for all 
initial conditions whenever E[log /(0, &)] = E[r$ — Pt] > 0. Alternatively, E[log /oo(£t)] = E[r t ] < 0 
implies asymptotic extinction with probability one for all initial conditions. Conditional persistence 
occurs when both of these inequalities are reversed. 

4. Positive and negative density dependence 

For populations exhibiting positive and negative density dependence, the fitness function f(x,£) 
can increase or decrease with density. For these general fitness functions, we prove several results 
about asymptotic extinction and persistence in the next two subsections. 

4.1. Extinction. We begin by showing that assumptions 
Al: Epog/(0,&)] <0, and 

A2: there exists 7 > 0 such that x 1— > f(x, £) is increasing on [0, 7) for all £ E E, 
implies asymptotic extinction occurs with positive probability for populations at low densities. 
Furthermore, we show this asymptotic extinction occurs with probability one for all positive initial 
conditions whenever the extinction set {0} is "accessible" i.e. there is always a positive probability 
of the population density getting arbitrarily small. More specifically, we say {0} is accessible from 
Ac [0, 00) if for any e > 0, there exists 7 > 0 such that 

F[3t > 0 : X t < e \ X 0 = x] > 7 
for all x £ A. We call a set B C M+ invariant if P [X x E B\X Q E B] = 1. 
Theorem 4.1. Assume Al and A2. Then for any 5 > 0, there exists e > 0 such that 

¥ [lim X t = 0\X 0 = x > I 

_t— >co 
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for all x G [0, e]. Furthermore, if {0} is accessible from [0, M) for some M > 0 (possibly +oo), and 
[0, M) is invariant, then 

P [lim X t = 0\X 0 = x] =1 

_t— >oo 

for all x G [0,M). 

There are two cases for which one can easily verify accessibility of {0}. First, suppose that 
f(x,£) = g(x)£. If (£t)i<o is a sequence of log-normal or gamma-distributed i.i.d. random variables 
and x i — y xg(x) is bounded (i.e. there exists M > 0 such that xg(x) < M for all x). Then, it 
follows immediately from the definition of accessibility that {0} is accessible from [0, oo). Hence, 
in this case E[log /(0, £t)] < 0 implies unconditional extinction. Since log-normal random variables 
and gamma random variables can take on any positive value, we view this case as the "large noise" 
scenario i.e. there is a positive probability of the log-population size changing by any amount. 

Alternatively, for sufficiently, small noise, there are a set of simple conditions for accessibility of 
{0}. Define F : R+ x E -)• R + by F(x,£) = f(x,£)x and the "unperturbed model" F 0 : R + -4 R + 
by Fq(x) = F(x,E[l;i\). For any x G R, define x + = max{0, x}. A system (1) satisfying the 
following hypotheses for e > 0 is an e -small noise system: 

HI: Fq is dissipative, i.e. there is a compact interval [0, M] and T > 1 such that Fq(x) G 

[0, M] for all x G [0, oo), 
H2 : F[F 0 (x) - £< F(x, 6) < F 0 (x) + e] = 1 for all x G M+, 

H3: for all x G R+ and all Borel sets U C [(^(a;) — e) + ,Fq(x) + e] with positive Lebesgue 
measure, there exist a > 0 and 7 > 0 such that 

P[F{z,Zi) G [/] >a 

for all z G [(x — "f) + ,x + 7]. 
The first assumption ensures that the unperturbed dynamics remain uniformly bounded. The 
second assumption implies that the noise is e-small, while the third assumption implies the noise 
is locally absolutely continuous. 

Proposition 4.2. Assume the difference equation xt+i = Fo(xt) has no positive attractor. Then 
there exists a decreasing function e : R+ — > 1R+ such that, for any M > 0, there exists an invariant 
set K D [0, M] such that {0} is accessible from K whenever (1) is an e{M)-small noise system. 

As a direct consequence of Theorem 4.1 and Proposition 4.2, we have 

Corollary 4.3. For any M > 0, there exists £0 > 0 such that if (1) is an e-small noise system for 
£ < £0; the dynamics induced by F$ has no positive attractor, and assumptions A 1-2 hold, then 

F [lim X t = 0|X 0 = 

_t— YOO 

for all x G [0,M]. 

4.2. Persistence. When E[log/(0, £1)] > 0 and there is only negative density dependence, The- 
orem 3.1 ensured the system is stochastically persistent. The following theorem shows that this 
criterion also is sufficient for models that account for negative and positive density dependence. 

Theorem 4.4. If 

(i) E[log/(0,6)] >0, and 

(ii) there exist x c > 0 such that E[sup| ;c>Xc | log/(x,£i)] < 0 and K[sup x<Xc log + f(x,£i)] < 00, 
then for all e > 0 there exists M > 0 such that 



limsupIT([l/M,M]) > 1 — e almost surely 

t— >oo 



whenever Xq = x > 0. 
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Remark 4.5. If there exists x c > 0 such that f(x,£) is a decreasing function in x on [x c ,oo), 
E[log f(x c , £i)] < 0, and K[sup x<x log + /(x, £i)] < oo ; then condition (it) in Theorem 4-1 is satis- 
fied. 

When the invasion criteria is not satisfied (i.e. E[log /(0, £i)] < 0), conditional persistence may 
still occur. For instance, suppose the stochastic dynamics have a positive invariant set B C (0, oo): 
there exists 7 > 0 such that B C [7,00) and F x [X t £ B for all t > 0] = 1 for all x € B. When 
such a positive invariant set exists, populations whose initial density lie in B persist. The following 
proposition implies that conditional persistence only occurs if there is such a positive invariant set. 

Proposition 4.6. Assume A1-A2. If (1) is bounded in [0,M] (i.e. F[X t < M for allt >Q] = l), 
then either lirm._j.oo Xt = 0 with probability one whenever Xq = x > 0, or there exists a positive 
invariant set B C (0, M] . 

In the case of small noise, the following proposition implies the existence of a positive attractor 
for the unperturbed dynamics is sufficient for the existence of a positive invariant set. In particular, 
conditional persistence is possible when E[log /(0, £1)] < 0. 

Proposition 4.7. Assume that A C (0, 00) is an attractor for the difference equation xt+i = Fo(xt). 
Then there exists a bounded positive invariant set K whenever (1) satisfies H2 for e > 0 sufficiently 
small. 

4.3. Mate-limitation and predator-saturation with negative density-dependence. To il- 
lustrate Theorems 4.1,4.4 and Propositions 4.2,4.7 we apply them to models accounting for negative 
density-dependence and positive density-dependence via mate-limitation or predator-saturation. 
The deterministic version of these models were analyzed by Schreiber [2003]. 

To account for negative density-dependence, we use a Ricker type equation. In the case of the 
mate-limitation model, the fitness function becomes 

x 

(2) f(x, £) = exp(r — a x )j^_ — where £ = (r, a, h) 

where r is the intrinsic rate of growth in the absence of mate-limitation, a measures the strength of 
infraspecific competition, and h is the half-saturation constant as described in Section 3.2. In the 
absence of stochastic variation in the parameters r, a, h, the dynamics of persistence and extinction 
come in three types [Schreiber, 2003]. If f(x,£) < 1 for all x > 0, then all initial conditions 
go asymptotically to extinction. If > 1 for some x > 0, then dynamics of extinction are 

governed by the smallest positive fixed point M and the critical point C of F(x) = xf(x,£). 
If F(F(C)) > M, then there is a positive attractor in the interval (M, 00) for the deterministic 
dynamics. Alternatively, if F(F(C)) < M, then the model exhibits essential extinction: asymptotic 
extinction occurs for Lebesgue almost every initial density, but there is an infinite number of 
unstable positive periodic orbits. In particular, there is no positive attractor. 

To account for environmental stochasticity, we assume, for illustrative purposes, that rt is uni- 
formly distributed on the interval [r — e, r + e] with r > 0 and 0 < e < r. Furthermore, we assume 
that a = 1 and h > 0. As E[log /(0, &)] = —00, Theorem 4.1 implies that limt_ 5 . OC) Xt = 0 with 
positive probability for initial conditions Xq sufficiently close to 0. When the deterministic dynam- 
ics support a positive attractor (i.e. F(F(C)) > M) and the noise is sufficiently small (i.e e > 0 
sufficiently small), Proposition 4.7 implies that the density Xt for the stochastic model remains in a 
positive compact interval contained in (M, 00). Alternatively, if the deterministic dynamics exhibit 
essential extinction and the noise is sufficiently small, Proposition 4.2 implies limt_>.oc Xt = 0 with 
probability one for all initial densities despite the deterministic dynamics having an infinite number 
of unstable periodic orbits. Finally, when e is sufficiently close to r (i.e. the noise is sufficiently 
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Figure 2. Asymptotic dynamics of extinction and persistence for the stochastic mate limitation 
model with negative density-dependence. For each parameter value, the model was simulated 
10, 000 time steps for multiple initial conditions. The final 1, 000 points of each simulation are 
plotted. Model details: The fitness function is f{x,£) = exp(r — ax)j^ where r is uniformly 

distributed on [4.5 — e, 4.5 + e] 



large), Theorem 4.1 implies that liim ; _ >00 Xt = 0 with probability one for all positive initial con- 
ditions. This later outcomes occurs whether or not the deterministic dynamics support a positive 
attractor. Each of these outcomes are illustrated in Fig. 2. 

For the predator-saturation model, we use the fitness function 

(3) f(x, £) = exp ( r — a x — ] where £ = (r, a, h, P) 



h + x 



where h and P are the half-saturation constant and the maximal predation rate, respectively, 
as described in Section 3.3. The dynamics of persistence and extinction for this model without 
stochastic variation come in four types [Schreiber, 2003]. If /(0,£) > 1, then there is a positive 
attractor whose basin contains all positive initial densities. If < 1 for all x > 0, then all 

initial conditions go asymptotically to extinction. If f(x, £) > 1 for some x > 0, then dynamics 
of extinction are governed by the smallest positive fixed point M and the critical point C of 
F{x) = xf(x,£). If F(F(C)) > M, then there is a positive attractor in the interval (M, oo) for 
the deterministic dynamics. Alternatively, if F(F(C)) < M, then the model exhibits essential 
extinction. 

To account for stochasticity, we assume for simplicity that Pt is uniformly distributed on the 
interval [P (1 — e),P(l + e)] for some P > 0 and 0 < e < 1. Furthermore, we assume that a = 1, 
r > 0, and h > 0. When E[log /(0, £f)] = r — P > 0, Theorem 4.4 implies the system is stochastically 
persistent. Alternatively, when E[log/(0,£t)] = r — P < 0, Theorem 4.1 implies that linn._-.oo Xt = 0 
with positive probability for initial conditions Xq sufficiently close to 0. Assume r < P. If the 
deterministic dynamics support a positive attractor (i.e. F(F(C)) > M) and the noise is sufficiently 
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Figure 3. Asymptotic dynamics of extinction and persistence for the stochastic predator 
saturation model with negative density-dependence. For each parameter value, the model was 
simulated 10,000 time steps for multiple initial conditions. The final 1,000 points of each 
simulation are plotted. Model details: fitness function /(_,£) — exp(4 — 4x — P t /(1/12 + a;)) with 

P uniformly distributed on P[l — e, 1 + e]. 



small (i.e. e > 0 sufficiently small), Proposition 4.7 implies that the density Xt for the stochastic 
model remains in a positive compact interval contained in (M, oo). Hence, the population exhibits 
conditional persistence. Alternatively, if the deterministic dynamics exhibit essential extinction 
and the noise is sufficiently small, Proposition 4.2 implies linn._-.oo Xt = 0 with probability one for 
all initial densities. Finally, when e is sufficiently close to 1 (i.e. the noise is sufficiently large) 
and P > r, Theorem 4.1 implies that Iim£_ s . 00 Xt = 0 with probability one for all positive initial 
conditions. Each of these outcomes are illustrated in Fig. 3. 



5. Discussion 

A demographic Allee effect occurs when individual fitness, at low densities, increases with popula- 
tion density. If individuals on average replace themselves at very low densities, then the population 
exhibits a weak Allee effect. Alternatively, if there is a critical density below which individuals do 
not replace themselves and above which where they do, then the population exhibits a strong Allee 
effect. It is frequently argued that environmental stochasticity coupled with a strong Allee effect can 
increase the likelihood of a population falling below the critical threshold, rendering them particu- 
larly vulnerable to extinction [Courchamp et al., 1999, Stephens et al., 1999]. While this conclusion 
is supported, in part, by mathematical and numerical analyses of stochastic differential equation 
models [Dennis, 2002, Krstic and Jovanovic, 2010, Yang and Jiang, 2011], these earlier analyses 
are specific to a modified Logistic growth model with Brownian fluctuations in the log population 
densities. Here, we analyzed discrete-time models allowing for general forms of density-dependent 
feedbacks and randomly fluctuating vital rates. Our analysis demonstrates that environmental 
stochasticity can convert weak Allee effects to strong Allee effects and that the risk of asymptotic 
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extinction with strong Allee effects depends on the interaction between density-dependent feedbacks 
and environmental stochasticity. 

When environmental fluctuations (£t) drive population dynamics (xt+i = f(xt,£t+i)%t), an Allee 
effect is best defined in terms of the geometric mean G(x) = exp(E[log f(x, ^)]) of fitness. If 
the geometric mean G(x) is an increasing function at low densities, an Allee effect occurs. If 
this geometric mean is greater than one at low densities (G(0) > 1), then we proved that the 
Allee effect is weak in that the population stochastically persists: the population densities spends 
arbitrarily little time at arbitrarily low densities. When the geometric mean is less than one at low 
densities (G(0) < 1), the stochastic Allee effect is strong: for populations starting at sufficiently low 
densities, the population density asymptotically approaches zero with positive probability. Since 
the geometric mean G(0) in general does not equal the intrinsic fitness /(0,E[^]) at the average 
environmental condition, environmental stochasticity can, in and of itself, shift weak Allee effects 
to strong Allee effects and vice versa. For example, a shift from a weak Allee effect to a strong 
Allee effect can occur when a population's predator has a fluctuating half-saturation constant. 
Specifically, for the predator saturation model considered here, the geometric mean at low densities 
equals G(Q) = r — E[P//it] where r is the intrinsic rate of growth of the focal population, P 
is proportional to the predator density, and h% is the fluctuating half saturation constant of the 
predator. As Jensen's inequality implies that G(0) < r — P/M[ht], fluctuations in h% can decrease 
the value of G(0) from > 1 to < 1 and thereby shift a weak Allee effect to a strong Allee effect. 

In the absence of negative density-dependent feedbacks, we proved that there is a dynamical 
trichotomy: asymptotic extinction for all initial densities, unbounded population growth for all 
positive initial conditions, or a strong Allee effect (i.e. G(0) < 1 but G(x) > 1 for sufficiently large 
x). When a strong Allee effect occurs and environmental fluctuations are large (i.e. the support of 
log f(x, £t) is the entire real line for all x > 0), populations either go asymptotically to extinction or 
grow without bound with probability one. Moreover, both outcomes occur with positive probability 
for all positive initial conditions. 

Liebhold and Bascompte [2003] used models with only positive density-dependence to examine 
numerically the joint effects of Allee effects, environmental stochasticity, and externally imposed 
mortality on the probability of successfully exterminating an invasive species. Their fitness function 
was 

/(*,&) = exp ( 7 (x-C) + 6) 

were C is the deterministic Allee threshold, 7 is the "intrinsic rate of natural increase," and £t are 
normal random variables with mean 0. Since G(0) = exp(— 7C) < 1 and lim x _ s . 00 G(x) = +00 for 
this model, our results imply both extinction and unbounded growth occur with positive probabil- 
ity and, thereby, provide a rigorous mathematical foundation for Liebhold and Bascompte [2003] 's 
numerical analysis. Consistent with our simulations of a stochastic mate limitation model, Liebhold 
and Bascompte [2003] 's found that the probability of persistence increases in a sigmoidal fashion 
with initial population density. In particular, environmental stochasticity increases the probabil- 
ity of persistence for populations initiated at low densities by pushing their densities above the 
deterministic Allee threshold. Conversely, for populations initiated at higher densities, environ- 
mental stochasticity can increase the risk of asymptotic extinction by pushing densities below this 
threshold. Indeed, we proved that the probability of asymptotic extinction approaches zero as 
initial population densities get large and the probability of asymptotic extinction approaches one 
as initial population densities get small. 

Since populations do not grow without bound, negative density-dependent feedbacks ultimately 
dominate population growth at higher population densities [Wolda and Dennis, 1993, Turchin, 
1995, Harrison and Cappuccino, 1995]. While stochastic persistence never occurs with a strong 
Allee effect, extinction need not occur with probability one. Whether or not extinction occurs 
for all positive initial densities with probability one depends on a delicate interplay between the 
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nonlinearities of the model and the form of environmental stochasticity. A sufficient condition for 
unconditional extinction (i.e. extinction with probability one for all initial conditions) is that the 
extinction set {0} is "attainable" from every population density state. Attainability roughly means 
that the population densities become arbitrarily small at some point in time with probability one. 
For populations whose densities remain bounded from above, we proved a dichotomy: either there 
exists a positive invariant set for the process or {0} is attainable in which case there is unconditional 
extinction. Whether this dichotomy extends to unbounded population state spaces remains an open 
problem. 

When environmental stochasticity is weak and there is a strong Allee effect, the "unperturbed" 
population dynamics determines whether extinction occurs for all initial conditions or not. By 
"weak" we mean that the unperturbed dynamics F are subject to small, compactly supported 
random perturbations (i.e. xt+\ — F(xt) lies in an interval [— e,e] for e > 0 small). The existence 
of a positive attractor is necessary for conditional persistence in the face of weak environmental 
stochasticity. This result confirms the consensus in the mathematical biology community, that the 
existence of a positive attractor ensures that population trajectories can remain bounded away 
form extinction in the presence of small perturbations [Schreiber, 2006]. 

For populations exhibiting a strong Allee effect and conditional persistence at low levels of envi- 
ronmental stochasticity, there is always a critical level of environmental stochasticity above which 
asymptotic extinction occurs with probability one for all initial population densities. Mathemati- 
cally, there is a transition from the extinction set {0} being inaccessible for part of the population 
state space at low levels of environmental stochasticity to {0} being accessible for the entire popu- 
lation state space at higher levels of environmental stochasticity. We have illustrated this transition 
in stochastic models of mate-limitation and predator-saturation with negative density-dependence. 
Surprisingly, for the predator-saturation models, our numerical results show that environmental 
stochasticity can lead to asymptotic extinction at intermediate predation rates despite conditional 
persistence occurring at higher and lower predation rates. This effect, most likely, is due to the 
opposing effects of predation on overcompensatory feedbacks and the Allee threshold resulting in 
a larger basin of attraction for the extinction state at intermediate predation rates. 

While our analysis provides some initial insights into the interactive effects of Allee effects and 
environmental stochasticity on asymptotic extinction risk, many challenges remain. Many popula- 
tions exhibit spatial, ontogenetic, social, or genetic structure. Proving multivariate analogs to the 
results proven here could provide insights on how population structure interacts with the effects 
considered here to determine population persistence or extinction. Furthermore, all populations 
consist of a finite number of individuals whose fates are partially uncorrelated. Hence, they ex- 
perience demographic as well as environmental stochasticity [Ackleh et al., 2007]. In accounting 
for bounded, finite population sizes in stochastic models, extinction in finite time is inevitable. 
However, these models often exhibit meta-stable behavior in which the populations persist for long 
periods of time despite both forms of stochasticity and Allee effects. This meta-stable behavior 
often is associated with quasi-stationary distributions of the finite-state models. Studying to what 
extent these distributions have well definite limits in an "infinite-population size" limit is likely to 
provide insights into these metastable behaviors [Faure and Schreiber, 2014] and provide a more 
rigorous framework to evaluate the joint effects of stochasticity and Allee effects on population 
persistence and ultimately their consequences for conservation and management. 
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about models with Allee effects and environmental stochasticity?" at a conference many years 
ago that inspired the work presented here. This research was supported in part by U.S. National 
Science Foundation Grants EF-0928987 and DMS-1022639 to Sebastian Schreiber. 



Downloaded from http://biorxiv.org/on September 18, 2014 



12 g. roth and s.j. schreiber 

6. Appendix 

This Appendix provides proofs of all the results in the main text. In Section 6.1, we prove a 
general convergence result based on an accessibility assumption that leads to the proofs of Theo- 
rems 4.1 in Section 6.2 and Theorem 3.2 in Section 6.3. In Section 6.4, we prove Proposition 4.2. 
In Section 6.5, we prove Proposition 4.6, and in Section 6.7, we prove Proposition 4.7. 

We begin with some useful definitions and notations. Let "B be the Borel a-algebra on M + . Let 
5,, denote a Dirac measure at i.e. 5y(A^ — 1 if y G A and 0 otherwise for any set ^4 G 23. Let 
R+ U {oo} be the one-point compactification of R+ and assume that {oo} is a fixed point for the 
system (1). For a sequence (xt)t>o C M+ U {oo}, we write xt > D when (xt)t>o converges to 

t— >oo 

D C M+ U {oo}, i.e. if for any neighborhood U of D, there exists T > 0 such that x% G U for all 
t > T. 

We consider the trajectory space formed by the product = IR+ equipped with the product a- 
algebra 23 N . For any x 6 M + (viewed as an initial condition of trajectory), there exists a probability 
measure F x on £1 satisfying 

F x [{u £n:uj 0 £A 0 ,...,uj k £ A k }\ = F[X 0 eA 0 ,...,X k € A k \X 0 = x] 

for any Borel sets Aq, . . . , A k C M+, and Px[{^ G : wo = x}\ = 1. The random variables Xt are 
the projection maps 

X t : n -> R+ 

OJ I y OJt- 

For the proof of Theorems 4.1 and 3.2, we consider the space of the environmental trajectories 
equipped with the product u-algebra £ N , and the probability measure Q on satisfying 

Q[{e G E N : e 0 G E 0 , . . . ,e k G E k }} = P[^ 0 G E 0 , . . . G ,B fc ] 

for any Borel sets ■ ■ • j E k C E. For now on, when we write e G we mean e = (et)t>o- Since 
E is a Polish space (i.e. separable completely metrizable topological space), the space E^ endowed 
with the product topology is Polish as well. Therefore, by the Kolmogorov consistency theorem, 
the probability measure Q is well defined. In this setting, the random variable £t is the projection 
map 

it: £ N R + 
e H> e t . 

We use the common notation E (resp. E x ) for the expectation with respect to the probability 
measure Q (resp. ¥ x ). 

Let x G M+ and Q x = {uj G f2 : ujq = x} be the cylinder of the trajectories starting at x. The 
continuous function tp : E N — > Q x defined component-wise by 

(4) p(e) t = X t ( v (e)) = xf(x, e 1 )f(X 1 (<p(e)),e 2 ) ■ ■ ■ /(^(^(e)), e t ), 

links the probability measures Q and ¥ x . In fact, the pushforward measure of Q by tp is the 
probability measure F x , i.e. for any Borel set A C Q x , F X (A) = Q(ip~ 1 (A)). 

Recall that a set A C M+ U {oo} is accessible from B C M+ U {oo} if for any neighborhood U of 
A, there exists 7 > 0 such that 

F x [3t>0 : X t eU}>~f 

for all x G B. A subset C C [0, 00) is invariant for the system (1) if P^-Xi G C] = 1 for all x G C, 
and it is positive if C C (0, 00). 
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6.1. Convergence result. 

Proposition 6.1. Let B C _i+ be an invariant subset for the system {!), and A C B be an 
accessible set from B. Assume that there exists 0 < 5 < 1 and a neighborhood U of A such that 



X t 



t— >OD 



-> A 



>l-6, 



for all x G U. Then 
for all x G B. 

Proof. Define the event C = {Xt 
U of A such that 



X t 



t— >oo 



-> A 



t— >oo 



> A}. By assumption there exists 5 > 0 and a neighborhood 



[e] > i - 5, 



for all x G U. Fix x G B and define the stopping time r = inf{i > 0 : Xt G {/}. Since A is 
accessible from 5, there exists 7 > 0 such that P x [t < 00] > 7. The strong Markov property 
implies that 



[e] > [ Px. [e] i {T 

Jn 

> (1-5)7, 



<OQ -jdP x 



The Levy zero-one law implies that linn._-.oo E^ [I g 1 3_ ] = le lP_-almost surely, where 3^ is the cr- 
algebra generated by {Xi, . . . , Xt}. On the other hand, the Markov property and invariance of B 
imply that E x [l e \J t ] = Px t [6] > (1 - 5)7- Hence F X [Q] = 1. 

□ 



6.2. Proof of Theorem 4.1. To prove the local extinction result, assume E[log/(0,^)] < 0 
and that there exists 7 > 0 such that x h-> /(_,£) is increasing on [0,7) for all £ £ E. Since 
E[log /(O, £1)] < 0 and x h-> f(x,£) is monotone on [0,7], there exists 0 < x* < 7 such that 
E[log/(x*,^i)] < 0 and /(_,£) < for all x G [0,s*) and all £ £ E. The Law of Large 

Numbers implies that 

1 <_1 

lim -V log /(_*,&) <0, 

t— >oo t £ — * 

s=0 

Q-almost surely. Define the random variable 

R = inf 1 



/Oo, .1) • ••/Ot-l, 
where the infimum is taken over the set Ut>i{( x o> • • • > #t-l) G [0, „*)*}. As 



(5) 



^ t-i t-i 
limsup - S~] log/ (a; s ,6+i) < lim log <0 

s=0 s=0 



Q-almost surely for any sequence {xt}t>o lying in [0, x*), -R > 0 Q-almost surely. Let F C £' N be the 
set of probability 1 for which the limits in (5) exist and R > 0. Choose e G T and x G [0,x*R(e)]. 
Let 93 : _v N — > Q x be the function defined by (4). The definition of R(e) implies by induction that 



Downloaded from http://biorxiv.org/on September 18, 2014 



14 G. ROTH AND S.J. SCHREIBER 

Xt(<p(e)) < 5 for all t > 1. Hence, our choice of x* implies that 

1 1 <_1 

limsup-logX t (< / 3(e)) = lim - Vlog(/(X s (^(e)), e s+1 ) 

+ — t t— >co t ' — ' 



t— >OD 



s=0 

t-l 



- } im 7y]log(/(^,e s +i)) 

t— >oo t ' — * 



s=0 



< 0 



for all e G T. As Q[T] = 1, 



11 1 

; [limsup - logX t = 0] = Qflimsup - log X t o ^ = 0] > Q[i? > — 

t— >oo I t— >oo t 



x*n 



for all n > 0 and x £ (0, jM. 

Fix 5 > 0. Since {i? > ^;}n>i is an increasing sequence of events and Q[U n {R > = 1, 

lim n _ >oc Q[i? > -4-] = 1 which implies that there exists N > 0 such that Q[R > -kw] > 1 — 5. 
Hence 









lim X t = 0 

.t— >oo 


> Q 






- X * N 



> 1-5, 

for all x G (0, 

Now assume {0} is accessible from [0, M) for some M > 0 (possibly +oo). Applying Proposi- 
tion 6.1 to A = {0} and -B = [0, M] (resp. -B = [0, oo)) implies lim^oo X t = 0 with probability one 
whenever Xq = x G [0, M). 

6.3. Proof of Theorem 3.2. To prove the extinction result, suppose that 7 = E[log /oo(£t)] < 0 
and fix x G M+. Let Let </? : £7 — >• be the function defined by (4). Since x h-> f(x,£) is an 
increasing function for all £ G £7, we have 

t-i t 
*tMe)) = Hf(X s (ip(e)),e s+1 )x < Y[foo(e s )x, 

s=0 s=l 

for all e G £^ N . By the Law of Large Numbers, 

limsupylogAiO(e)) < lim - ( V] log /oo(e s ) + logx ) = 7<0 

for Q-almost all e G i? N . Therefore 

P x .[lim X t = 0] = Q[lim X t o p = 0] = 1, 

t— »oo t— >oo 

which completes the proof of the first assertion. 

To prove the unbounded growth result, suppose that a = E[log /(0, £t)] > 0 and fix x G (0, 00). 
Since x h-» f(x,£) is an increasing function for all ( 6 £, we have 

f-i f 
*t(¥>(e)) = II f(X 8 (<p(e)), e s+1 )x > [] /(0, e,)x 



s=0 



s=l 



for all e G E^. By the Law of Large Numbers, 



lim inf - log X t (ip(e)) > lim - ( V) log /(0, e a ) + log x ) 



a > 0 



Downloaded from http://biorxiv.org/on September 18, 2014 



ALLEE EFFECTS AND ENVIRONMENTAL STOCHASTICITY 

-almost all e E E^. Therefore 



15 



5 x [Iim X t = oo] 

t— >oo 



[ lim Xtotp 
t— >oo 



OO 



] = 1, 



which completes the proof of the first assertion. 

To prove the Allee effect result, fix 0 < 5 < 1 and assume that E[log/(0, £i)] < 0 and 
E[log /oo(£i)] > 0. The first assertion of Theorem 4.1 implies that there exists m > 0 such that 



(6) 



lim X t = 0 

_t— >oo 



> 1-5, 



for all x E (0,m]. To prove the second part of the result, consider the process Y t = ^ 
to the event {Xq > 0}. It satisfies the following stochastic difference equation: 

Y t+1 =g(Y t ^ t+l ), 



conditioned 



where g : M+ 

9M = tAb 



X Q — > K is the continuous function defined by g(y,£) 
. By definition of Yt, we have, for all m > 0, 



for all y > 0 and 



lim Y t = 0 



t— >OD 



lim X t 

t— >oo 



OO 



for all y E (0, m]. 

Since, E[log<?(0,£i)] = — E[log /oo(£i)] < 0, the first assertion of Theorem 4.1 applied to {Yt)t>0 
implies that for any 0 < 5 < 1, there exist L > 0 such that 



lim Yj 

t— >oo 



0 



> 1-5, 



for all y E [0, L). Therefore, for any 0 < 5 < 1, there exists M = 4- > 0 such that 



(7) 



lim Xt = cxd 

t— >oo 



> 1-5, 



for all x E [M, oo). 

To prove the last part, assume that {0,oo} is accessible from 
there exist m, M > 0 such that 



and fix 5 > 0. By (6) and (7), 



X t ► {0, oo} 



t— >oo 



> 1-5, 



for all a; E [0, m] U [M, oo). Proposition 6.1 applied to A = {0, oo}, B 
concludes the proof. 



and U = [0,m]U[M,oo) 



6.4. Proof of Proposition 4.2. The proof consists of combining two deterministic arguments 
with a probabilistic argument. The three of them use the concept of (e,T)-chain introduced by 
Conley [1978]. An (e,T)-chain from x to y in M, for a mapping Fq : IR+ — > M+, is a sequence of 
points xq = x, xi, . . . , xt-i = y in IR + such that for any s = 0, . . . ,T — 2, — Fo(xj)| < e. x 

chains to y if for any e > 0 and T > 2 there exists an (e, T)-chain from x to y. 

The following Propositions are the deterministic ingredients of the proof and are proved in 
Schreiber [2006]. 

Proposition 6.2. Let A be an attractor with basin of attraction 23(A) and V C U be neighborhoods 
of A such that the closure U of U is compact and contained in 23(A). Then there exists T > 0 and 
5 > 0 such that every 5 chain of length t >T starting in U ends in V . 

Proposition 6.3. If Fq : — > IR+ satisfies HI and has no positive attractor, then for all x E M+, 

e > 0 and T > 0 there exists an e chain from x to 0 of length at least T. 

The probabilistic ingredient is an adaptation of Proposition 3 in Schreiber [2007] to our frame- 
work. 
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Proposition 6.4. Assume (1) satisfies H3 for eq > 0. If x G M + chains to 0, then, for all e < Eq, 

there exists a neighborhood U e of x and f3 £ > 0 such that 

F z [3t > 0 : X t < e] > /3 e 

for all z G U £ . 

Proof. For any a, b G M + , define I(a, b) := [(a — 5) + , a + b]. Let e < £q and xq = x, x\, . . . , xt = 0 
be an |-chain from x to 0. There exists 7t > 0 such that It := I(Fo(xt-i),'jt) C [0, e]. Assumption 
H3 implies that there exist 7t_i > 0 and at > 0 such that := I(a;t-ij7t-i) C J(i*b(aJt-2)> e) 
and 

P^LYi G J t ] > at 

for all z G it-i. Since Jt_i C I(Fo(xt-2), e), assumption H3 implies that there exist 74-2 > 0 and 
at-i > 0 such that 

P Z [X! G 7 t -l] > at-l 

for all z G /*_2 := I(xt-2,~ft-2)- Repeating this argument, there exist Io,...,It C M+ and 
ai, . . . , at > 0 such that for all s = 1, . . . , t 

Vz[Xi G I,] > a s 

for all z G / s _ 1. Define a := mm s a s . The Markov property implies that, for all z G M+, 

P z [Xt G It] = E z [ . . . E Xt _ 2 [P*^ [Xi G It] 
Since for all s = 1, . . . , t, F Z [X± G I s ] > ati a _ 1 , 

F z [X t G > a 1 for all z£l 0 . 
Choosing U = Iq and j3 = a competes the proof of the proposition. □ 

Lemma 6.5. Let e > 0 and V <ZU be bounded subsets o/M+. Assume that the system (1) satisfies 
H2 /or e and i/iai i/iere exists T > 0 suc/t f/taf every e chain of length t >T starting in U ends in 
V. Then there exists a bounded invariant set K D U for the system (1). Moreover, if U C (0, 00), 
then K C (0,oo). 

Proof. Assume that U C (0, 00). Let L = sup xg (y{Fo(x)} + svcp x&u {x] + Te. Since every e chain 
of length t > T starting in U ends in V, assumption H2 implies that, for any x G U, Xt G (0, L) 
for all i > 0 with probability one whenever Xq = x. Define the positive bounded Borel set 

K := {x G (0, L) : F x [3t > 0 : X t > L] = 0}. 

Since U C K, K is nonempty. To show that -fC is invariant for (1), let x G K. By the Markov 
property, 

0 = F x [3t > 0 : X t > L] > E x F y [3t > 0 : X t > L]l K c 

Since, F y [3t > 0 : X t > L] > 0 for all y G /sT c , F x [Xi G ^sT c ] = 0. Hence, K is a positive invariant 
set for (1). 

If 0 G U, then it follows from the same arguments that {x G [0, L) : P a [3t > 0 : Xt > L] = 
0} D C/ is invariant for (1). □ 

Proof of Proposition 4.2. Since i*o is dissipative, there exists an attractor A such that 
2(A) = R + . Let V C M+ be a neighborhood of A and M 0 = inf{M > 0 : V C [0,M]}. For 
any M > Mq, Proposition 6.2 applies to A, V and [0, M]. Hence there exists e : [Afo,oo) — > R+ a 
decreasing function, T : [Mo, 00) — > N such that, for every M > Mq, every e(M) chain of length 
t > T(M) starting in [0, M] ends in V. We extend the functions e to IR + by defining e(M) = e(Mo) 
for all M < M 0 . 
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Fix M > 0, and assume that (1) is an e(M)-small noise system. If M > M$, then Lemma 6.5 
implies that there exists an invariant set Km D [0, M] for the system (1). 

Assume that Fq has no positive attractor. Propositions 6.3 and 6.4 imply that, for all x G Km 
(the closure of Km) and all e < e(M), there exists a neighborhood U Xj£ of x and (3 Xt£ > 0 such that 

¥ z [3t >0 : X t <e] >/3 Xj£ 

for all 2 G U x>£ . Compactness of Km implies that, for any e < e(M), there is /3 e > 0 such that 

¥ z [3t > 0 : X t < e] > /3 e 

for all z G -f^M- Hence {0} is accessible from -Km- 

If M < M 0 , then [0, M] C [0, M 0 ] C K Mo - Moreover, by definition of e(M), i^ A / 0 is invariant for 
(1) and {0} is accessible from Km 0 - This concludes the proof. 

6.5. Proof of Theorem 4.4. After showing that the system (1) is almost surely bounded, the 
almost surely persistence follows as in the proof of Theorem 1 in Schreiber et al. [2011]. Define 

V : M+ -> M+ by V(x) = x, a, /3 : Q -¥ R+ by a(f) = sup{/(x,£) : x > x c } and /3(f) = 
sup{xf(x,£) : x < x c }. Hence, for any f G $7 and x G M+, we have V(xf(x,£)) < a(£)V(x) + /?(£)• 
By assumption, E[lna] < 0 and E[ln + /3] < oo where ln + (x) = max{0,x}. Hence Proposition 4 in 
Benaim and Schreiber [2009] implies that the system (1) is almost surely bounded. 

6.6. Proof of Proposition 4.6. Assume E[log /(0, £i)] < 0. Recall, we say (1) is unconditionally 
extinct if lirn^oo Xt = 0 with probability one whenever Xq = x > 0. If there exists a positive 
invariant set B C (0, M], then the system (1) can not be unconditionally extinct. If the system (1) 
is not unconditionally extinct, then, by Theorem 4.1, {0} is not accessible from [0, M]. Therefore 
there exists e > 0 such that for all n > 1, there exists x n G [0, M\ such that 

¥ Xn [3t >0X t <e}<-. 

n 

Since the event {3t > 0 X t < e} is an open set of £1, compactness of [0, M] and weak* continuity 
of x I— > Pa, imply there exists x G [0, M] such that P^pi > 0 Xt < e] = 0. Define the non empty 
positive set B = {x G [0,M] : F x [3t > 0 X t < e] = 0} C [e,M]. To show that B is invariant, let 
x G B. We will show that P^f-Xi G B] = 1. By the Markov property, 

0 = F x [3t >0X t <e]> E x [PxJ^ > 0 s.t. X t < e]t B 4X 1 

Since P y [3t > 0 s.t. X t < e] > 0 for all y G B c , P x [Xi e B c ] = 0 for all x e B. Hence, £ is 
invariant. 

6.7. Proof of Proposition 4.7. Assume that A C (0,oo) is a positive attractor with basin of 
attraction B{A). Let V C U C B(A) be positive compact neighborhoods of A. Proposition 6.2 
applies to A, V and U. Hence, there exists e and T > 0 such that every e chain of length t > T 
starting in U ends in V. Assume the system (1) satisfies H.2 for e. Then Lemma 6.5 implies that 
there exists a positive bounded invariant set for (1) which concludes the proof. 
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